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ABSTRACT 

The observed Fe n(UV+optical)/Mg n AA2796,2804 flux ratio from a gravitationally lensed quasar 
B1422+231 at z=3.62 is interpreted in terms of detailed modeling of photoionization and chemical 
enrichment in the broad-line region (BLR) of the host galaxy. The delayed iron enrichment by Type la 
supernovae is used as a cosmic clock. Our standard model, which matches the Fe n/Mg n ratio, requires 
the age of 1.5 Gyr for B1422+231 with a lower bound of 1.3 Gyr, which exceeds the expansion age of 
the Einstcin-de Sitter f2o = 1 universe at a redshift of 3.62 for any value of the Hubble constant in the 
currently accepted range, i?o=60-80 kms -1 Mpc -1 . This problem of an age discrepancy at z = 3.62 
can be unraveled in a low-density £Iq ^ 0.2 universe, either with or without a cosmological constant, 
depending on the allowable redshift range of galaxy formation. However, whether the cosmological 
constant is a required option in modern cosmology awaits a thorough understanding of line transfer 
processes in the BLRs. 

Subject headings: cosmology: theory — galaxies: evolution — quasars: emission lines — quasars: 
individual (B1422+231) 



Mg ii AA2796,2804 doublet lines, so that the ratio of their 
fluxes, emergent from the same ionized zone, is a direct 
indicator of the abundance ratio (Wills, Netzer, & Wills 
1985). These wavelengths are redshifted to the near in- 
frared at z 3, and quasars are the only objects that are 
bright enough for their spectra to be studied. In this re- 
gard, Kawara et al.'s (1996) infrared observation of the 
Fe ii to Mg ii flux ratio for the gravitationally lensed 
quasar B1422+231 at z = 3.62 is particularly important. 
In this Letter, we compare their data with detailed mod- 
els of chemical enrichment, and attempt to determine the 
age of the broad-line region (BLR) of B1422+231, which 
allows us to set a combined constraint on the cosmological 
parameters and the epoch of galaxy formation. 

2. BASIC EQUATIONS FOR CHEMICAL EVOLUTION 

The strength of Fe n emission from quasars at various 
redshifts implies a much more significant iron supply from 
SNe la to the gas in the BLR of quasars as compared to the 
solar neighborhood (Wills et al. 1985; Elston, Thompson, 
& Hill 1994; Kawara et al. 1996). Such an iron-rich gas 
results if quasar host galaxies are associated with an initial 
burst of star formation and then the ejected iron from SNe 
la is restored to the gas after the burst and the cessation of 
star formation. This situation is explored by modeling the 
chemical enrichment with two basic ingredients (Tinsley 
1980) such as (1) the star formation rate (SFR) propor- 
tional to some power of the gas fraction C(t) = Vk[f g (t)] k , 
and (2) the initial stellar mass function (IMF) having a 
time-invariant mass spectrum <p(jn)dm oc m~ x dm normal- 
ized to unity between the lower and upper mass limits 
(m h m u ). 

The SFR coefficient for quasar host galaxies has been 
constrained in a narrow range from a variety of their ob- 
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1. INTRODUCTION 

Ages of distant objects at high redshift, z > 1, unam- 
biguously constrain the expansion age of the universe at 
that redshift, and can be used to determine the ultimate 
fate of the universe (Turner 1991; Kennicutt 1996). There 
is a growing interest in use of Type la and II supernovae 
(SNe la and II) as a nucleosynthesis-clock in interpreting 
the abundance pattern of heavy elements seen in the spec- 
tra of high-redshift galaxies (Hamman & Ferland 1993; 
Matteucci & Padovani 1993). 

Extensive calculations of explosive nucleosynthesis show 
that SNe II are the sites of alpha and iron production, 
whereas SNe la produce a large amount of iron relative to 
alpha elements (Nomoto, Thielemann, & Wheeler 1984). 
Since progenitors of SNe la have a lifetime of t\ a ~ 1 Gyr, 
longer by two or three orders of magnitude than SNe II 
(Truran 1987), the switchover of the iron source from SNe 
II to SNe la creates a break in the alpha/iron abundance 
ratio at a time at which a significant number of SNe la 
start to explode. Accordingly, the nucleosynthesis clock 
records such a switchover and assigns an age above or be- 
low ii a ~ 1 Gyr, depending on whether the alpha/iron 
ratio is smaller than or equal to that of SN II origin, re- 
spectively. This feature acts like an alarm clock set at 
ii a ~ 1 Gyr, and motivates us to measure the alpha/iron 
ratio at z £ 3, to obtain a firm constraint on the expan- 
sion age of the universe to be compared with that derive 
from the Hubble constant, for which we adopt Hq =60-80 
km s -1 Mpc -1 (Freedman 1998) as the current best esti- 
mate. 

Among various alpha elements, magnesium is ideal 
for comparing with iron, because the best-studied range 
of rest-frame wavelengths from 1500 to 6000A includes 
the emission features of Fe 11 multiplet lines and 
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served features, that is, fk=i — 7.5 Gyr 1 from their 
metal abundances (Padovani & Matteucci 1993; Matteucci 
& Padovani 1993) and v k=1 = 6.7 - 7.6 Gyr" 1 from the 
N v/C iv and N v/He n line ratios of high-redshift quasars 
(Hamann & Ferland 1992, 1993) having high mctallicities 
like B1422+231. We, therefore, adopt the higher value 
of Vk = \ = 7.6 Gyr" 1 as the standard value in this pa- 
per. Theoretical arguments indicate that the IMF origi- 
nates from fragmentation of the gas cloud which occurs 
almost independently of local physics in the gas (Low 
& Lynden-Bell 1976- Silk 1977). ' A solar-neighborhood 
IMF would therefore be a good approximation, and we 
adopt a Salpeter slope of x = 1.35 and a mass range from 
mi = 0.05 to m u = 5OM (Tsujimoto et al. 1997). Then 
the gas fraction f g (t) and the heavy-element abundance 
Zi(t) in the gas change with time according to 
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respectively, where nit is the turnoff mass when the main- 
sequence lifetime, t mi is equal to time t, r(m) the fraction 
of the ejected material from a star of mass m, r w (m) the 
fraction of the ejected material without newly synthesized 
elements from that star, and yi the heavy-element yield 
from an SN II or la. Since all these stellar quantities, ei- 
ther calibrated or constrained by nearby stars, should also 
apply to any galaxy, we can use the formula of Renzini 
& Buzzoni (1986) for m t and the updated nucleosynthesis 
calculations by Nomoto's group (Tsujimoto et al. 1995) 
for high-redshift galaxies. 

The fraction of stars that eventually produce SNe la is 
A = 0.055 for 3 — 8M Q and ^4 = outside this mass range. 
The lifetime of their progenitors is t\ a w 1.5 Gyr and its 
possible spread is modeled using the power-law distribu- 
tion function g(ii a ) <x tj a (7 > 0) which is bounded in a 
specified range of t\ a and normalized to unity. These ba- 
sic quantities A and t\ a for SNe la have been constrained 
mainly from the observed [O/Fe] break at [Fe/H]^ —1 
in the solar neighborhood (Yoshii, Tsujimoto, & Nomoto 
1996). From theoretical considerations Kobayashi et al. 
(1998) recently proposed that SNe la should occur only 
when the progenitor's metallicity is above a critical value 
of [Fe/H]~ -1. While this implies that the [O/Fe] break 
at [Fe/H]^ —1 may only weakly constrain t\ ai the intro- 
duction of a critical metallicity for SNe la resolves an ob- 
servational puzzle (Carney 1996): The Galactic globular 
clusters (GCs), spanning an age difference of several Gyrs, 
exhibit a constant alpha/iron abundance ratio of the SN 
II origin all the way from [Fe/H]^ —2 to —0.5 without 

4 The synthetic spectrum with a A _1 -typc extinction law fitted to the data gives a dust extinction below Ay 
that this effect hardly changes our measurement of the Fc Il/Mg II flux ratio. 



a break. In fact, with this critical metallicity, the GC 
data place even firmer limits of t\ a = 1 — 3 Gyrs from the 
box-shaped (/-distribution with 7 = (Tsujimoto & Yoshii 
1998). We therefore use this constraint on t\ a to date the 
delayed enrichment of iron relative to alpha elements in 
high-redshift galaxies. 

3. AGE DATING FOR THE B1422+23I SYSTEM 

The rest-frame UV-optical spectrum of B1422+231 was 
obtained with the IR detector array mounted on the 
KPNO 4m telescope by Kawara et al. (1996). The ob- 
served Fe n(UV+optical)/Mg 11 AA2796,2804 flux ratio is 
12.2 ± 3.9, 4 which is comparable to 8.9 derived from 
the composite spectrum of intermediate-redshift quasars 
at z = 1 — 2 (the Large Bright Quasar Survey (LBQS), 
Francis et al. 1991) and 7.8 ± 2.6 for nine low-redshift 
quasars at z = 0.15 — 0.63 (Wills et al. 1985). Using pho- 
toionization models where important excitation processes 
are all included, Wills et al. (1985) showed that within the 
wide range of physical parameters in their solar-abundance 
models the Fe 11/Mg 11 flux ratio is firmly bounded be- 
tween 1.5 and 4 with a typical value of 3 for the likely 
optical depth r(Bac)= 0.5 — 1 of the hydrogen Balmer 
continuum. Therefore, the observed strength of Fe 11 emis- 
sion from their sample yields an overabundance of Fe by 
a factor of 3 with respect to Mg (see also Netzer & Wills 
1983; Collin-Souffrin, Hameury, & Joly 1988; Graham). 
Accordingly, with a reasonable assumption of Fe 11/Mg 11 
oc Fe/Mg, the flux ratio for B1422+231 leads to an 
abundance ratio of [Mg/Fe]= -0.61^;^ ([Mg/Fc] = 
log(Mg/Fe)— log(Mg/Fe) Q ), where the first errors quoted 
correspond to the uncertainties in transforming the flux to 
abundance ratio and the second to the observational errors 
in measuring the flux ratio. Throughout this paper we use 
the transformation formula [/(Fe n)/J(Mg 11)]© = 31 1 5 . 
We note that its uncertainty, especially in the direction of 
giving a larger [Mg/Fe] (or smaller age) of B1422+231 is 
easily exceeded by the combined uncertainties of the flux 
measurement and chemical evolution model which will be 
discussed below (see Table 1). 

Figure 1 shows the Fe 11/Mg 11 flux ratio (upper panel) 
and the logarithmic iron abundance [Fe/H] (lower panel) 
as a function of time in units of Gyrs. Using the Salpeter 
IMF and the nucleosynthesis prescriptions as described 
above, we calculate the standard model (vk=i = 7.6 
Gyr -1 ; thick line) together with the lower SFR model 
{vk=i = 6.7 Gyr -1 ; dashed line). The Fe 11/Mg 11 flux ra- 
tio is initially maintained at a low level reflecting the low 
Fe/Mg abundance ratio in the SN II ejecta, and then starts 
to increase with the enhanced Fe supply due to the onset of 
SNe la from 1 Gyr until 3 Gyr. Thereafter, the Fe 11/Mg 11 
flux ratio declines, because the metal-deficient gas having 
the low Fe/Mg abundance ratio of genuine SN II origin is 
released in the interstellar matter from the surface of low- 
mass stars, with lifetimes greater than 3 Gyr, that were 
born in the initial burst of star formation. The horizontal 
line in the upper panel is the observed Fe 11/Mg 11 flux 
ratio for B1422+231, and the shaded region brackets the 
range of the errors. 

The intersection between the theoretical curve by the 

0.1 mag for Bf 422+231, so 
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standard model and the horizontal line representing the 
observed Fe n/Mg n flux ratio in the upper panel gives 
the age of 1.50 Gyr for B1422+231, which, according the 
lower panel, requires [Fe/H]=+0.96 at that age. This 
high metallicity of ten times solar agrees well with the 
result derived by Hamann & Ferland (1992, 1993) from 
the N v/C iv and N v/He n line ratios of high-redshift 
QSOs with z = 2 - 4. 

We have examined the uncertainty of the estimated age 
and metallicity by repeating the calculations with differ- 
ent values of input parameters from our standard choice. 
Table 1 tabulates At = t(chd)— t(std) for each of the 
changed parameters and shows that the standard model 
already uses the parameter values which result in a rela- 
tively low estimate of the age. Use of a higher A(SN la) 
in the model makes this age smaller by ~ 0.1 Gyr, thus 
imposing a lower bound of 1.3 Gyr, if the uncertainties of 
the observed Fe n/Mg n flux ratio is taken into account. 
This result can be applied unless A(SN la) in the BLRs is 
drastically larger than the value constrained in the solar 
neighborhood. 

4. COMPARISON WITH THE EXPANSION AGE OF THE 
UNIVERSE 

Consider a galaxy at redshift z, which is assumed to 
form at zp. The age of this galaxy is the time taken by 
the universe to expand from zp to z. Because of its high 
redshift of z = 3.62, the lower age bound of 1.3 Gyr for 
B1422+231 constrains not only the cosmological param- 
eters but also the epoch of galaxy formation in a highly 
correlated manner. Figure 2 shows the (Qo,Ho) region 
allowed by the existence of B1422+231. Two values of 
zp — 7 and 10 are assumed for either a universe with 
Aq = (thin lines) or a flat universe with ilo + Ao = 1 
(thick lines). Evidently, the Einstcin-de Sitter universe 
(flo = 1, Ao = 0) requires H < 40 kms -1 Mpc -1 and 
zp £ 10. However, the currently accepted range for Hq 
does not allow a value as small as 40 kms -1 Mpc -1 for 

Figure 3 shows the predicted age-redshift relations for 
three representative combinations of the cosmological pa- 
rameters: the Einstcin-dc Sitter universe (f2o = 1, Ao = 0; 
left panel), an open universe (f2 = 0.2, A = 0; mid- 
dle panel), and a flat, A-dominated universe (Qo = 0.2, 
A = 0.8; right panel). Solid lines in each panel denote 
the result for Hq= 60 (upper line) and 80 (lower line) 
kms -1 Mpc , assuming zp — 10. For the purpose of com- 
parison, the result for zp = 5 and Hq= 60 kms -1 Mpc -1 
is shown by the dashed line. 

Figures 2 and 3 indicate that the age discrepancy, now 
encountered at z = 3.62 in the Einstein-de Sitter universe, 
is alleviated either in an open, Slo ~ 0.2 universe with 
zp £ 10, or in a flat, f2 ~ 0.2 universe with zp £ 6 — 7. 
With an age set at the lower bound of 1.3 Gyr, we obtain 
two solutions for two different ranges of zp. However, with 
an age of 1.5 Gyr as obtained with our standard model, 
an open universe becomes difficult unless fi ~ 0.1 and 
z F > 10. 

We note that the IMF in the star-bursting phase is con- 
sidered to have a shallower slope (Contini, Davoust, & 
Considere 1995) and/or a larger mi (Rieke et al. 1993) 
than in the quiescent solar neighborhood. Such an IMF 



weighted towards massive stars was also derived from the 
N/C ratio in the BLRs (Hamman & Ferland 1992, 1993). 
If this were the case in B1422+231, a longer age, exceeding 
even 1.5 Gyr, is required to reach the observed Fe n/Mg u 
flux ratio (see Table 1). The existence of such an old galaxy 
at z = 3.62 would rule out an open universe, and nonzero- 
A would have to be invoked. More Fe n/Mg n observations 
in high-redshift QSOs are evidently needed. 

The interesting possibility of a requirement for nonzcro- 
A has recently received some support from new distance 
determinations to SNe la discovered near z — 1 (Pcrlmut- 
ter et al. 1998; Garnavich et al. 1998). 

5. CONCLUSION 

Detailed models of chemical enrichment are used to esti- 
mate the age of a gravitationally lensed quasar B1422+231 
(z = 3.62) from the Fe n/Mg n flux ratio. Our standard 
model, which matches the data, requires an age of 1.5 
Gyr with a lower bound of 1.3 Gyr, whereas the Einstein- 
de Sitter universe (f?o = 1, A = 0) is 0.6-0.8 Gyr old 
at that redshift with current estimates of the Hubble con- 
stant Ho — 60 — 80 kms -1 Mpc -1 . This problem of the age 
discrepancy, now encountered at z = 3.62, was first iden- 
tified at z w by dating the GCs (15.8±2.1 Gyr, Bolte 
& Hogan 1995; 14.56±2.49 Gyr, Chaboyer et al. 1996) 
and the radioactive element 232 Th in a very old halo star 
(CS22892-052) (15.2 ± 3.7 Gyr, Cowan et al. 1997), and 
then at z = 1.55 by dating the UV-optical spectral en- 
ergy distribution for a high-redshift radio galaxy LBDS 
53W091 ( > 3.5 Gyr, Spinrad et al. 1997). 

In finding potential ways to remove the age discrep- 
ancy of the GCs, a brighter RR Lyare calibration of 
My(RR)« +0.2 mag, and therefore younger GC ages, 
was inferred from Hipparcos measurements of parallaxes of 
Galactic Cepheids (Feast & Catchpole 1997) and Galactic 
metal-poor subdwarfs (Reid 1997). However, an analysis 
of Hipparcos proper motions and parallaxes of Galactic RR 
Lyraes by Tsujimoto, Miyamoto, & Yoshii (1998) directly 
confirmed My(RR)« +0.6 mag, as previously accepted. 
Therefore, the age discrepancy of the GCs is not yet re- 
solved. 

The age discrepancy at z — 3.62 in this paper is based 
on the Fe/Mg abundance ratio derived from a simple scal- 
ing of their emission line fluxes and therefore needs to 
be confirmed by more detailed modeling of the emission 
line transfer of Fe n in the BLR using the newly available 
atomic data (e.g., Sigut & Pradhan 1998). In this regard, 
although the present analysis already suggests a nonzero-A 
universe (see Fig. 3), a definitive discrimination between 
an open universe and a nonzero-A universe awaits future 
efforts of elaborating the photoionization model as well 
as observing the Fe n/Mg n flux ratio in more QSOs at 
z ~ 3 — 5, because no other age-dating methods except for 
the nucleosynthesis-clock are able to be employed at such 
a great distance and early epoch in the universe. 
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Tabic 1 

Sensitivity of the estimated age and metallicity to the change of the model parameters 

t(std) = 1.50 Gyr [Fc/HJ(std)= +0.96 



Changes At A [Fe/Hj 

[7(Fc)//(Mg)J Q = 3 -> 1.5/4 +1.05/ - 0.13 +0.32/ - 0.10 

SFR v k=1 = 7.6 -> 6.7 +0.08 +0.01 

SFR k = 1 -> 2 a +0.63 -0.01 

IMF x = 1.35 1.1 +1.15 +0.39 

IMF to; = 0.05 -> 0.1 +0.26 +0.08 

SN la A = 0.055 -> 0.08 -0.14 +0.04 

SN la ff-form: 7 = -> l b +0.21 -0.04 



The coefficient Vk=2 is chosen so that the gas fraction at 0.5 Gyr should be the same as that from Vk=i = 7.6 Gyr" 
b The power-law index 7 = 1 corresponds to a case that SN la progenitors of longer lifetime, which are possibly less 
massive, are given a higher weight like the IMF. 
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Fig. 1. — Evolutionary behaviors of the Fcll/Mgll flux ratio (upper panel) and the logarithmic iron abundance [Fe/H] (lower panel) for 
theoretical models calculated with different SFRs. Shown arc the standard model with v k — j = 7.6 Gyr -1 (thick line) and the lower SFR 
model with Vk=i = 6.7 Gyr (dashed line). The horizontal line is the value of the Fcll/Mgll flux ratio measured for B1422+231, and the 
shaded area is the uncertainty in this measurement. 
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Fig. 2. — Constraints on the cosmological parameters (Ho, Qq, Ao) and the epoch of galaxy formation zp. Shown are the Qq versus Ho 
relations derived from the lower age bound of 1.3 Gyr for B1422+231 at z = 3.62. Thin lines represent the cases with zp = 7 and 10 for a 
universe with Ao = 0, and thick lines are those for a fiat universe with Qo + Ao = 1. The area below each line is the allowed reigion where the 
expansion age of the universe exceeds the lower age bound of 1.3 Gyr. Note that current estimates of the Hubble constant are in the range 
of H = 60 - 80 kms" 1 Mpc" 1 . 
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Fig. 3. — The age-redshift relations predicted for the Einstein-de Sitter universe or a flat, matter-dominated universe with (f!o, Ao) = (1, 0) 
(left panel), an open universe with (f2o,Ao) = (0.2,0) (middle panel), and a flat, A-dominatcd universe with (fio,Ao) = (0.8,0.2) (right 
panel). In each panel shown are the predictions for zp = 10 with Hq= 60 and 80 kms -1 Mpc -1 (solid lines) and for zp = 5 and Hq= 60 
kms _1 Mpc _1 (dashed line). The age of 1.5 Gyr for B1422+231 estimated by the standard model of chemical enrichment is indicated at 
z = 3.62 by the filled circle together with its lower age bound of 1.3 Gyr. As a summary of the existing age determinations at z ss 0, the 
ages of GCs (15.8±2.1 Gyr, Bolte & Hogan 1995; 14.56±2.49 Gyr, Chaboyer et al. 1996) are shown by the open circles and the age of the 
CS22892-052 star (15.2 ± 3.7 Gyr, Cowan ct al. 1997) by the cross. The age of LBDS 53W091 ( £ 3.5 Gyr, Spinrad ct al. 1997) is also 
indicated at z = 1.55. 



